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The study started in Ref. [l6l ] about the Dilaton mean field stabilization thanks to the effective 
potential generated by the existence of massive fermions, is here extended. Three loop corrections 
are evaluated in addition to the previously calculated two loop terms. The results indicate that the 
Dilaton vacuum field tend to be fixed at a high value close to the Planck scale, in accordance with 
the need for predicting Einstein gravity from string theory. The mass of the Dilaton is evaluated to 
be also a high value close to the Planck mass, which implies the absence of Dilaton scalar signals 
in modern cosmological observations. These properties arise when the fermion mass is chosen to 
be either at a lower bound corresponding to the top quark mass, or alternatively, at a very much 
higher value assumed to be in the grand unification energy range. One of the three 3-loop terms 
is exactly evaluated in terms of Master integrals. The other two graphs are however evaluated 
in their leading logarithm correction in the perturbative expansion. The calculation of the non 
leading logarithmic contribution and the inclusion of higher loops terms could made more precise 
the numerical estimates of the vacuum field value and masses, but seemingly are expected not 
to change the qualitative behavior obtained. The validity of the here employed Yukawa model 
approximation is argued for small value of the fermion masses with respect to the Planck one. A 
correction to the two loop calculation done in the previous work is here underlined. 
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I. INTRODUCTION 



The Dilaton is an essential ingredient of superstring theory, and constitutes a scalar field partner of the graviton 
Therefore, the background fields associated with the vacuum state of superstring theory should involve this 
field in common with the metric in the basic action. This is referred to as Dilaton gravity 0, Q. To the lowest 
level of approximation the Dilaton is a free and massless scalar field with a special kind of coupling to the matter 
fields. As a consequence of this coupling, a time varying Dilaton field determines time-dependent coupling constants. 
In order to overcome this difficulty the Dilaton should remain constant during the present stage of evolution of the 
Universe. Moreover, unless the Dilaton turns out to be very massive, its existence could lead to an observable "Fifth 
force" similar to the ones which are currently associated to the observations of the Dark Matter. The constraints 
posed by current experimental observations determine the lower bound on the mass of the Dilaton to be of the order 
to < KT 12 GeV H (but see @ for an attempt to make a running Dilaton consistent with late time cosmology). 

The Dilaton stabilization problem has been at the center of an intense research activity in recent times because 
of its physical relevance. It should be emphasized that the Dilaton is one of various scalar fields appearing in the 
formulation of superstring theory in the low-energy limit. The sizes and shapes of the extra spatial dimensions 
associated with superstring theory are also leading to additional scalar fields, called "moduli fields" . The stabilization 
of such moduli fields has been the object of recent attention particularly in connection with Type IIB superstring 
theory. The introduction of fluxes within the compactification spaces has made it possible to stabilize various moduli 
fields 0- Also, gaugino condensation Q has been employed to stabilize the Dilaton field in the context of heterotic 
superstring theory Q and in string gas cosmology 10]. 

It should be remarked that, since Dilaton stabilization has special relevance for late time cosmology, there is 
motivation for finding mechanisms which do not directly rest on the concrete assumptions defining the nature of the 
extra dimensions. An additional motivation to search for alternative Dilaton stabilization mechanisms comes from 
String Gas Cosmology (SGC). The SGC [THHl] is a model of early universe cosmology which employs new degrees of 
freedom and symmetries of string theory, and couples these elements with gravity and Dilaton fields into a classical 
action background model. The Universe is considered to start as a compact space containing a gas of strings. Since in 
string theory there is a maximal temperature for a gas of closed strings, the initial state of the cosmological evolution 
in SGC will be a phase of almost constant temperature, the so called "Hagedorn phase". The SGC is able to define 
a non-singular cosmology in which there is no starting Big Bang explosion. It has been noted that the thermal 



2 



fluctuations in a gas of closed st ring s in the Hagedorn phase can justify the scale-invariant spectrum of cosmolo gica l 
fluctuations observed in Nature [l3l [l4j , with a particular prediction of a slight blue tilt for gravitational waves [l5j j . 
However, the consistency of the picture requires that the Dilaton field be fixed during the Hagedorn phase. Therefore, 
in the SGC theory the Dilaton needs to be fixed at very early times and at very late times. 

Thus, clarifying the mechanisms of Dilaton field stabilization is an important question in particle physics today. It 
is worth noting that the universal type of coupling of the Dilaton to the matter fields not only leads to an unwanted 
effect as the time-dependence of the coupling constants but it also furnishes the possibility that quantum effects due 
to the interaction of the Dilaton with matter might generate interesting contributions to the effective potential of the 
Dilaton. 

In a previous work published in Ref. (l6j , we started to explore this question. The work considered the cosmological 
periods when the additional spatial dimensions of superstring theory were already stabilized and the study was done 
in the framework of a four-dimensional field theory. The objective of study was then the interaction of the Dilaton 
with massive fermions. Such masses can be defined by fluxes about internal manifolds. In late time cosmology, the 
masses could had been generated after supersymmetry breaking. In an alternative early universe cosmology, one 
may consider thermally generated fermion masses. In Ref. [16[ it was considered a simple form for the Dilaton 
gravity action in which a massive Dirac fermion term was added 17[. The action was chosen in the Einstein frame, 
which does not show any Dilaton field dependence in the kinetic terms for the fermions. On the other hand, the 
fermion mass becomes a function of the Dilaton, involving a universal exponential factor in Dilaton gravity @, H| ■ 
The chosen action described the low energy effective interaction of Super- Yang-Mills fermions with the Dilaton field 
in superstring theory [l6[ . The effective potential for the Dilaton field was evaluated up to two loop corrections in the 
small Dilaton radiative quantum field limit. That leads to a Yukawa like interaction term which allows standard QFT 
calculations. A fixed value of the cosmological scale factor was assumed. The outcome of the work was, thanks to the 
appearing of logarithms in the loop calculations, that the Dilaton field appeared in the result in powers multiplied by 
the exponential factors of the field. This structure, in the one loop approximation clearly indicated the spontaneous 
generation of vacuum mean value of the Dilaton field. 

Motivated by the dynamical generation of the Dilaton result in Ref. [lq . we here will address the evaluation of 
next corrections 3-loop terms to the 2-loop evaluation of the effective potential for the Dilaton field. The main issue 
to be explored is the possibility of the appearance in the improved potential of the stabilizing effect which were in 
fact absent in the second order correction, and which are suspected to be created by the existence of massive matter 
upon the mean value of the Dilaton. 

The results obtained, at least indicate, for the fermion mass being selected at the GUT or the top quark mass 
scales, that the mean value of the Dilaton field tends to be stabilized at a high value being close to the Planck mass 
or the GUT scale, respectively. Therefore, it is suggested that the appearance of mass for matter in the course of 
the evolution of the Universe can generate a stabilizing action on the vacuum expectation value of the Dilaton field 
making it unobservable. This effect will tend to stop the time evolution of the mean value, as it is convenient for 
String Theory consistency. 

The work also present an study of the validity of the linear approximation of the Dilaton exponential factor which 
leads to the simpler effective Yukawa theory employed. It follows that for the two values of the large fermion masses 
assumed (the top quark and GUT scale ones) the approximation should work well, after assuming that the low energy 
effective effective action of string theory is the bare one in the renormalization of the model. A brief resume about 
the procedure employed to estimate the Dilaton mass is also given. 

We can point out, that in the process extending the work to include higher loop corrections, we have noticed that 
in Ref. [lj| the kinetic term of the Dilaton Lagrangian was chosen with a negative sign. This selection, although not 
changing the one loop correction, led to a sign change of the 2-loop terms, which suggested the existence of minima 
in the effective action argued in Ref. (l6| . However, in spite of this non physical adopted assumption in that work, 
the indication about the dynamical generation in Ref.flU] remained a valid one, because the change in the metric did 
not affected the one- loop correction, the basic quantity indicating the dynamical generation effect. The present work 
corrects the result for the two loop terms, and indicates that its place in the stabilizing effect over the Dilaton field 
is played by higher order contributions. 

The paper proceeds as follows: In Section II, the notation and basic formulation are given. Section III presents the 
elements of the three loops evaluation of the effective potential. Section IV discuss the results of the calculation. 
In Section V the results are resumed and commented. Finally, Appendix A presents the investigation of the Yukawa 
model approximation and the review on the scheme for evaluating the generated Dilaton mass. 
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II. THE DILATON ACTION AND GENERATING FUNCTIONAL 



Let us consider a model of the Dilaton field interacting with fermion matter in the form 



S = J d 4 x^/^gjx)(^g^(x)d ft <j> r (x)d v <l> r (x) + $(a;)(i Qg 9 v - m)*(x) 

-V(x)g Y 4> T \x)^{x) + j{x)4> r {x) + ^(x)rj(x) +r](x)^(x)), 
m = exp(a* <j))mf, 



g Y = a m, 
3 



a = — - 



(x°,x\x 2 ,x 3 ), d =d-d, {l^,l v } = 2g llv {x), 
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That is, we are considering the Dilaton field interacting with a massive fermion in the Einstein frame, in which the 
metric g^ v has been approximated by the Minkowski metric in order to simplify the evaluation. The gravitational 
constant is here explicitly introduced, and natural units are employed for the distances and mass. The vacuum value 
of the Dilaton field is named as (f> and its radiative part is called (jf . Note that we are assuming the radiative part is 
small in order to retain only the first term in the expansion of the exponential. This is the Yukawa approximation 
which is here employed. In appendix A it is argued that it can be a a good approximation for the two values of the 
large fermion masses considered here : the top quark mass and a GUT scale one, assumed that the the ratio between 
the fermion mass and the Planck one is very much smaller than one. All the results will be functions of the vacuum 
field <j) an d the fermion mass m. 

The parameter defining the Dilaton field dependent exponential, the Planck length n — lp and mass mp are defined 
by the expressions 



8irGti 



1 



G 
h 

c 



l P = — = 8.10009 x 10 

mp 

6.67 x 10~ 8 cm 3 g" 1 s~ 2 , 
1.05457 x 10~ 27 cm 2 g s" 1 , 
2.9979245800 x 10 10 cm s" 1 



cm, 



(7) 

(8) 

(9) 
(10) 
(11) 



In the above formula for the action, the coordinates and times are measured in cm, the masses m in the natural 
unit cm -1 and the Dilaton field is dimensionless. 

Starting from the classical action, we will consider a 3-loop correction to the effective action, assuming a homogenous 
and time independent value of the Dilaton mean field <\> as 



U(4) 



= -v en (0), 



(12) 



where is the four dimensional volume. In order to eliminate the explicit appearance of the gravitational constant 
from the diagram technique for evaluating the effective action, we could absorb it by redefining the Dilaton field value 
and the a* constant as 
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<P/k, 



a = a k = ~~^ K i 



9y 



g Y K. 



(13) 
(14) 
(15) 



After these changes, the above written classical action S, to be used for generating the Feynman expansion can be 
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expressed as follows 

S^,*,<p r ,<p] = J d 4 x (±g^(x)d»<p r (x)d v <p r (x) + *(x)(i 9 ^ 9 v - - mMx) 

-H{x)g Y (p r (x)y(x) +j(x)( ip + ip r (x)) + \F(ar) r)(x) +rj(x) *(ar)). (16) 

The expansion is considered in d = 4 — 2e dimensions for implementing dimensional regularization scheme. Ac- 
cordingly, the coupling constant g Y should be modified by the introduction of the regularization scale parameter \x 
as follows 

9v = ^{9vf, 
where g Y is the usual coupling constant in four dimensions. 



A. The generating functional and the effective action 



In this subsection, for the sake of definiteness, we will sketch the main expressions defining the perturbative 
calculation to be considered in what follows. The generating functional of the Green functions Z , its connected 
part W and the mean field values will defined by the formulae 



z [v,vJ] 

W[rj,r],j] 

S W 
i 5j(x) 

5 W 
i 5r](x) 
S W 
—i 6r)(x) 



V^V^Vy r exp(i S [*, tf , <p r , if}), 
^log Z[rj,r),j], 
<p+ (ip r (x)), 

<*(*)>, 
(*(*)>■ 



(17) 
(18) 
(19) 
(20) 
(21) 



Note that the mean Dilaton field tp is considered as homogeneous and the mean value of the radiative part (ip r (x)} 
will be assumed to vanish when the sources are zero. The effective action is defined as the Legendre transform of Z 
depending on the mean field values as: 



r[<*),<*), v + <^)] 
s r 

s r 

ST 



-Aog Z[r}, V ,j} - J dx\j(x){<p + (<p r (x))) + n(x) +rj(x) (*(x)) ], (22) 

-3(x), (23) 

- V (x), (24) 

r}(x). (25) 



The expression for Z, after writing the Yukawa vertex part of the Lagrangian in terms of the functional derivatives 
over the sources and integrating the gaussian functional integral that remains, leads to the Wick expansion formula: 



z [v,V,j] = exp[i J dx 9Y~j: 



exp 



D(x-y) = 



Sj(x) —iSr](x) iSr](x) 
J dx dy (rj(x)S(x - y)r](y) + ^j(x)D(x - y)j{y)) 
dp d cxp(— i p.(x — yj) 



) d m- ■yVpn 
dk d exp(— i k.(x — y)) 
(2ir) d -(fc 2 -ie) : 



(26) 
(27) 
(28) 
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in which S and D are the fermion and Dilaton free propagators, respectively. The notation for fermions and scalar 
field related quantities, and the definition of the Feynman rules for the generation of the analytic expressions for the 
various contributions are exactly the ones described in Ref. (l8| |. for the cases of scalar and fermion fields. Specifically, 
for the momentum space rules, the propagators and the only existing vertex are graphically illustrated in figure [T] 



r • >- 



4> Hnr) 



p 1+ p 2 



FIG. 1: The figure illustrates the Feynman rules for the particular Yukawa model approximation adopted for the Dilaton 
action 



III. EFFECTIVE POTENTIAL EVALUATION 



Let us present in this section the evaluations of the effective potential for the Dilaton field following after employing 
the perturbative expansion described in the past section. The diagrams to be considered are depicted in Fig. [2] 
They include up two three loop corrections. The contributions will be exactly evaluated for the one and two loops. 
In addition, the three loop term D32 also can be analytically calculated in terms of Master integrals. However, the 

three loop diagrams D31 and D 33 are here determined only in their leading terms of order log (j^j . We expect to 

be able in evaluating the non leading corrections (lower powers of log (jfj) in extending the work. Let us discuss 
the results for each diagram in various subsections below. 




FIG. 2: The one, two and three loops Feynman diagrams considered in this work. The one and two loop corrections Di and 
D2 are exactly calculated. In the case of the three loops terms, the D32 is completely evaluated in terms of the listed Master 
integrals in Ref. [TjJ. The D31 and D33 were determined only in their leading logarithm correction. 
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A. One loop term Di 

The analytic expression for the one loop diagram D\ and its derivative over m 2 have the forms 

= J ^-Trlog(m 2 - P 2 ), (29) 
" rF« = «*0 / tS^^- (30) 



d m 2 J (2ir) d i m 2 — p 

The result for the momentum integral entering in the derivative of I^ 1 ) over m 2 , after divided by fj, 2e (in order 
to define a 4-dimensional energy density) and integrated over m 2 , allows to write for the one loop effective action 
density (See Ref. pjj) 

r (1) A/ m, 9e 87T 2 ~ e . 

= J^FW = m y (2^p* r(_1 + £) - (31) 

After employing the minimal substraction (MS) scheme, that is, getting the finite part by eliminating the pure pole 
part in e the Laurent expansion of 7(771, e) and taking the limit e — > 0, the one loop contribution to the effective 
action density as a function of m and fi becomes 



71 (to, h) = 0.0506606m 4 f 2. log ( ^ J - 2.95381^) . (32) 

Note that the negative of this term, which gives the one loop effective potential leads to a the dynamical generation 
of the Dilaton field for positive values of a* <fi as follows from log(m) = log(m,f) +a* cj>. This was the effect which 
motivated the study started in Ref. [16]. 

B. Two loop term D2 

For the two loop contribution D 2 the analytic expression is 

T (2) _1, , 2 f dpf dpi 4{m 2 + Pl .p 2 ) 



, \ T( ' 1,0,2 f d Pl d t 



) d i (m 2 — p\)(m 2 — pf)(Pl — P2) 2 



1 („o \22d-4 [ dqf dq d 4(l + gi.«a) 

-{g Y j to 



2 VJY ' J (2*)*i(2*)*i(l-qi)(l 

dqf dq d 1 



2(.g° ) 2 TO 4 m" 4e (2 



(2n)*i (2n) d i (1 - q{)(l - q 2 ){ qi - q 2 f 



\d^i-.7T^n 03) 



2 K J (2^)Ml-g 2 )- 

where the identity q\ .q 2 = \{q 2 — 1 + q 2 2 ~ 1) + 1 — 5(91 — 92) 2 has been used. The two momentum integrals appearing 
in the last line are the simplest Master integrals for scalar fields as listed in Ref. [l9(. The results for them in that 
reference are: 

dqf dq d 2 1 _ (d-2)(n) d T(l-l) 2 



(2n) d t (2n) d i (1 - q 2 ){\ - q 2 )( qi - q 2 ) 2 2(d - 3)(2n) 2d ' 

dqf 1 W'r(i-|) 



(2n) d i {l-q 2 ) (2Tr) d 
They allow to write for the regularized two loop effective action density the expression 



(35) 



4/ m N Af 2(o° ) 2 M d , d-2 h / d\ 2 , . 
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Expanding in Laurent series in e and disregarding the pole part in the limit e — > 0, leads to the two loop 
perturbative contribution to the effective action 

y 2 (m,n) = 0.0000200507(5?- ) 2 m 4 (48.log 2 ( — ) - 173.783 log— + 183.83 ). (37) 

As it was noticed in the Introduction, in Ref. [16| it was employed an inappropriate negative kinetic term for the 
Dilaton field. This change, although not affecting the one fermion loop contribution, which is not altered by the sign 
of the boson propagator, drastically modified the sign of the two loop term which linearly depends on the Dilaton 
propagator. In the previous evaluation, the two loop terms determined the existence of minima for the Dilaton 
potential. Therefore, the consequence of the change in sign fixed by the here consideration of the correct positive 
kinetic energy term, should be investigated in connection with the existence of stabilizing minima for the scalar field. 
This circumstance determined the motivation for the new three loop corrections considered in this work. 



C. Three loops terms 



Let us now consider the three loop terms. 



D 32 

The -D32 term is the only of the 3-loops diagrams which is not composed of two fermion or boson self energy 
insertions connected in series. For the D31 and D33 cases we had difficulties in reducing their contributions to a 
linear combination of tabulated Master integrals. This obstacle only allowed us to calculate their leading term in the 
expansion in log(^). However, for D 32 it was possible to express it as a sum over the Master integrals given in Ref. 
[l9j . The analytic expression of the diagram is 

r ( 32 ) = _ V (d)i {gY ^ f dpj dpj dpj 



4 y ' J (2ir) d i (2ir) d i (2ir) d i 
Tr[(m + p^7 M )(m + (p% + p% - p^ )7 M )(m +P3^ fJ ,)(m + Pi7 M )] 



(to 2 -p{)(m 2 - pj)(m 2 ~pl)(m 2 - (p 2 + pz - pi) 2 ){pi -P3) 2 {pi - Pi) 2 

_V(«0 1(^)4 f d P 2 



4 y ' J (2ir) d i (2?r) d i (27r) rf i 

to 4 + c 1 (p 1 ,p 2 ,P3)m 2 + c 2 (pi,p 2 ,P3) , 38 , 

(to 2 -p\)(m 2 -pl)(m 2 -pl){m 2 - (p 2 + p 3 ~ Pi) 2 )(pi -P3) 2 {pi - Pi) 2 ' 

ci(pi,Pi,P3) = 3pi-P3+Pi-Pi+Pi-P3+pI+pI~pI (39) 

Ci{pi,Pi,P3) = P? Pi -P3 + p\ Pi -P3 + p\ Pi -Pi - 2 Pi -Pi Pi -P3 ■ (40) 



Defining now 



2 2 
zi = p 1 -m , 

z 2 = p\- to 2 , 

Z3 = pl~ TO 2 , 



= [pi - Pi) , 

Z5 = (pi -P3) 2 , 

z% = {pi ~Pi+ P3) 2 - m 2 , (41) 

and employing various vectorial identities expressing the squares of the differences between any two momenta in terms 
of the scalar product between them and the squares of the considered momenta, the integral defining r^ 32 ) can be 
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written as follows 



r (32) = _y{d)h gY) A 



dpi d P d 2 dpi 

x 



A y "' ' J (27r) d i (27r) d i (27r) d z 
m 4 + ci(z) m 2 + C2(z) 

Zl Z 2 Z 3 Zi Z 5 ZQ 

z = [zi,z 2 ,Z 3 , 24, z 5 ,z 6 ), (42) 
3 

c i( z ) = 2 ( z i + Z2+ z 3 + ^6)-2(z 4 +z 5 ) +6m 2 , (43) 

c 2 (z) = -(ziZ6 + z 2 ^3 - z 4 z 5 + m 2 (zi + z 2 + z 3 + z 6 ) + 2m 4 ). (44) 

Therefore, there are one or two z factors in the denominator that can be canceled by the terms of the quadratic 
polynomial in these quantities. This allows the integral to be decomposed in a linear combination of the Master 
integrals listed in Ref. [19]. The result for action density 

p(32) 

732 K^e) = jpyffi, 
can be expressed in terms of only five of them as follows 

7 32 (m, M ,e) = -( 3 o ) 4 ™ 4 (^) 6 (W) + 8/ 2 ( e )-4I 3 (e) + J 5 (e)-:^ 
where the functions h(e), /2(e), hi 6 ), H e ) and 17(e) result to be given by 



h{e) 



2 _3(4-2 e )-9 7r -|(4-2 e)(5(4 _ 2e) _ l 8 ) Ml ( e )3 

l-2e + 

2 -3(4-2 e )- 67r -3(4-2 e ) ( 3 (4 _ 2e) _ 10)(3(4 _ 2e) _ 8) ^ Mg(e) „ 2(4 _8 2 e e) _ 7 M 4 (e)) 



(45) 



2 -3(4-2e)-2 -3(4-26) / M ( \3 (9 _ Of) 2 \ 

He) = j-^ ( j 2e + (3(4 " 2£) " 8)M4(£) ) ' (46) 

2- 3 ( 4 - 2£ )- 3 ^- 3 ( 4 - 2£ ) /2(2-2e) 2 M 1 (e) 3 x , .A 

He) = M x _ 2e + (3(4 - 2e) - 8)M 5 ( e )J , (47) 

4(e) = (2^)- 3 ( 4 - 2e )M 4 (e), (48) 

I 7 (e) = (2iT)- 3 ^M 5 (e), (49) 

in terms of the Master integrals (See Ref. [Hj]): 

Mi(e) = ^ 4 - 2e >rQ(2 e -4) + l) , (50) 

M 2 (e) = J*-^ M f)\ (51) 

v y 2(1 -2e) ' v 7 

M 3 (e) = 2^ 2 - 4 )r Q(4 - 2e)) T Q(2e - 1)) M 4 ( e ) 2 , (52) 

2 1 - 2e r(i(8-3(4-2e)))r(i(2e-l)) , _ 

M 4 (e = ttj— ^- 2 — ^ Z^.- V2V r-^-M^e 3 , (53 

W r(i(7-2(4-2e)))r(|(2e-2)) u; ' 1 ' 

M 5 (e) = (-2-|e - ie 2 + ^e 3 + 1(163 - 128C(3))e 4 + 

( + HPT + ^ " l0 S( 2 ) 2 )( 321 °g( 2 ) 2 ) " 168 C(3) (54) 
-256Li 4 (i) )e 5 ) Af 4 (e) 3 , 
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where the special functions Li n (i) and £(n) are defined as 

Lin (a:) = ( 55 ) 



2 k k r ' 
fe=i 



oo 

C(») = E^- ^ 



k=i fc 

Finally, the application of the before described MS procedure leads to the following formula for the contribution 
to the vacuum effective action density of the diagram D 32 

/ 777 \ / I YY1 \ 1 771 

,0 \4„„4 ii\ nnnnqom 1 a i^J> I l nnnnmcnri/iiA^ 4 I i n rainrnficssi ' 



732 (m, /i) = (ffy) 4 m 4 (0.00003291 14 log s ^— J - 0.000105904 log 4 (^— J + 0.0000165851 log J ^— J + 

0.000441159 log 2 - 0.00074347 log [^j + 0.000388237). (57) 

It can be noted that this term has a high quintic power of log 5 (^) which is also determined by the high pole of the 
e expansion present in the function I\. This is the highest power of the log (jfj expansion appearing in the results. 
The next higher power, the fourth one, also is arising in this term. 

2. Diagram Z) 31 

We were not able to exactly evaluate this contribution (and also the one associated to D 33 ) in terms of Master 
integrals. Therefore, for both of these terms we here limited ourself to evaluate their leading terms in the expansion 

in powers of log (jjj-^j- F° r this purpose, the use was made of the circumstance that (at variance with -D32, but in 

coincidence with Z?3 3 ) this term corresponds to a loop formed by two one loop self-energy insertions. Since these 
self-energy terms are explicitly calculable in terms of hypergeometric functions, both terms can be expressed as single 
momentum integral in d dimensions. The diagram has the original analytic expression 



r(3 i) = _ V (d)h av) 4. f dpj dpj dpj 

2 KyY > J (2n) d i (2n) d i (2i:) d 



(2it) d i (2it) d i (2w) d i 
Tr [(to + P2Jn)(m + PiJf_c)(m + p^^){m + p^)] 
(to 2 -p\) 2 (m 2 -p\)(m 2 -pf)(pi -p 3 ) 2 (pi -P2) 2 



dpi dpi dpi 



2^ ' J (2ir) d i (2ir) d i (2n) d i 
to 4 + d 1 (p l7 p 2 ,p 3 )m 2 + d 2 (pi,p 2 ,P3) 
(to 2 -p 2 ) 2 (to 2 -p 2 )(m 2 ~pf){pi -p 3 ) 2 (pi -P2) 2 ' 
dl(pi,P2,P3) = p\ + 2pi.p 2 + 2pi.p 3 +P2-P3 , (59) 
d2(Pl,P2,P3) = 2 pi.p 2 P1-P3 - P\P2-P3- (60) 



We define now the fermion self-energy integral and its related vector as follows 
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s 3 i(p 2 



d P i 



{2ir) d i (m 2 — p\)(pi — p) 2 
~ (£^ r(£) Jo dX X " £(m2 " p2(1 X) ~ lSY 



„ 2-Pi(l - e,e,2- e, - 
' r( £ )(m 2 )- — Li^, (61) 



(27r) d wv ' e-l 

/ 2n /* d Pl Pip 

W J 7 (27r)^(m 2 -p2)( pi _p)2 

= a(p 2 ) P„, (62) 

/ 2x P 2 +m 2 2 L(m,e) 
a(p ) = o S3i(p ) , (63) 



im ' e = Trhrr~2 — 2T= r i ~o • 64 

J (2n) a i (m z — p\) (2TT) a 2 

In the above expressions, the Feynman parametric integral was explicitly evaluated by employing the algebraic 
calculation program Mathematica. After performing the Wick rotation in the external and integration momenta and 
extracting the d-dimensional solid angle arising form the angular integrals, the expression for the action density can 
be written as follows 



p(31) 



get the expression 



Kg° Y ) 4 (fy 6£ c(m,e) 
e 2 

f(p, e) = e 2 h{p, e)r(e) 2 + f 2 (p, e) e T(e) + f 3 (p, e), (66) 



73i(m,/i,e) = / — —^f(p,e)dp, (65) 

6 Jo [p z + 1) 



/ 2 (p,c) = (2-^— ^) S * 1 (p 2 , £ )i*( £ ), (68) 



a_z 2 )V* ^ 

2p 2 
4p2 



/3(P.e) = - ( 69 ) 



*5i(P,e) = —J (70) 

c(m,e) = tt;— -r , (71) 

V ' T(i(4-2e)) V ' 

L*(e) = eL(l,e). (72) 

As it was mentioned before, we were not able yet to find an epsilon expansion (rigorous or sufficiently approximated 
numerical one) allowing to exactly evaluate this integral after removing the regularization. Therefore, in order to 
determine an approximation for 731 we have made use of an assumption suggested by an exploration done about the 
asymptotic power expansion at infinity of the integrand as a function of the momentum integration variable p. It 
followed that all the terms of the expansion after integrated, show a single pole structure in their Laurent expansion 
in e. Then, it suggests that the full divergence of the integral at d = 4 is defined by a single pole in e. Assuming 
this property, the extraction of the leading correction in log(^) should be defined by the maximal power of log(^) 
appearing in the coefficient of the zero order term in the expansion of the modified integral 

-6 



>(<#) 4 (^) c(m,0) 



oc p3 - 2e 



73i(m,/i,e) = / — — —^f(p,0)dp. (73) 

e Jo (p z + 1) 
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Note that any other power of e in the expansions of c(to, e) and f(p, e) will reduce the maximal order of the negative 
powers of epsilon in the full expansion of 731 (to, /i, e) which determines the leading correction in the expansion. For 
f(p, 0) it follows 



p 4 I7p 2 ; 7 1 

1024tt 4 ~ 1024tt 4 + 256tt 4 ~ 256ttV 

Then, the use of the formula 



/(P,0) = (74) 



I 



00 3-2£+m ^ n 

*- ndp = - ~(m-2e + 2)csc(~(m-2e)), (75) 

/o (p 2 + l) 2 4 V y V 2 V v y 

which shows the ^ singularity, allows to write for 731 the leading logarithm correction to its finite part 

731 (m,/i) = -0.0000228551( 5 ^) 4 to 4 log 3 (^) . (76) 



3. Diagram D33 

As it was remarked, this terms will be treated in a similar way as it was done for -D31. Now, the corresponding 
self-energy insertions will be the boson ones. Again, the two self-energy loops are explicitly calculable in terms of 
hypergeomctric functions. The starting analytic expression of the diagram is 



d P d dpf dp d 2 1 



4 var; J (2ir) d i (2ir) d i (2ir) d i (p 2 ) 2 
Tr[(m + p^7 M )(m + (p + pi)' y 7 i ,)]Tr[(m + p^l^jm + (gg + pYlv)] 
(m 2 — p 2 ) 2 (m 2 — (pi + p) 2 )(m 2 — p^) 2 {m 2 — (p% + p) 2 ) 



(27r) d i (27r) d i (27r) d i (p 2 ) 2 

(to 2 +pi-(pi + p))(to 2 +P2-{P2 +p)) 

(to 2 — p 2 ) 2 (to 2 — (pi + p) 2 ){m 2 — p 2 ,) 2 ^ 2 — (jp2 + p) 2 ) ' 



(77) 



where the fermion traces were evaluated for writing the second form of the integral. The last expression evidences 
the decomposition in two serial self-energy terms. 

After rotating to Euclidean space the momenta variables of the integration regions and the external momentum, 
the fermion selfenergy integral and its related vector integral can be written as follows (See Ref. [l8| ) 



•S33(<7 2 ,e) = 



dqf 1 



(2ir) d (m 2 + q1)(m 2 + (q + gi) 2 

^&(e) f'dxil + ^xil-x))- 
(4tt)2 Jo m 

(47T)* TO 2 
1 

,2„ 



T{q 2 ) = / dx (1 + q 2 x(l - x))- 



2— + y^T4)(l ^ ^fc)' 2 Fi(l - e, £ , 2 - e, j^fc + 1)) 

g(e-l) 



V33 M 2 ) = 



E - x (9 - vV + 4)(l + v4=t) £ ^1(1 - e, 6, 2 - e, §(—7*- + 1)) 
y<r+4 y<r+4 

g(e-l) 

jgi Pip 

(27r) d i (to 2 — p 2 )(m 2 — (p + pi) 2 



a(p 2 ) = ~^S33(p 2 ,e)- 
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Again the result for parametric Feynman integral was analytically evaluated thanks to the use of the algebraic 
calculation program Mathematica. 

Thus, after extracting the Euclidean angular integrals and performing some transformations, the density of the 
effective action 

r (33) 

7 33 (m,^,e) 



can be expressed as a single momentum integral in the range (0,oo) as follows 



-6e 



4c(m,e)(g Y ) 4 (fX ,<*, p3 -2e 

7 33 (m,Ai,e) = / dp — —^g{p,e), (78) 

e Jo (p z + r z j 

g{p,e) = £ 2 9l (p,#) 2 lj2(p,f)fr( £ ) + j3 (p,e), (79) 

9i{m,e) = (y + 2 ) 4 3 (P,ef, (80) 

g 2 (m, e) = -2 (y + 2^ L(e) S * 3 (p, e), (81) 

g 3 (m,e) = L(e) 2 , (82) 

SjhM = r ^ . (83) 

c(m,e) = 77- -t , (84) 

r(i(4-2e)) 

i(£) = (2^ r( " 1 + £) - (85) 

Finally, by employing a similar procedure for extracting the leading logarithmic correction in log (jjf^j f° r D 3 \, the 
analogous contribution for Z? 33 follows in the form 

733 (m, M ) = -0.000329114 (g Y ) 4 m 4 log 3 (— I . 



IV. DISCUSSION 



Lets us now comment the results obtained in previous sections for the effective action density. The total effective 
potential value v(m, /i), is given by the sum of all the evaluated terms after changing their sign. The total potential 
and its various contributions are written below 



v(m, fi) 


= vi(m, ju) + \ 


v\{m,n) 


7i(to,/z) 


TO 4 


m 4 


v 2 (m, /i) 


72(m,/z) 


m 4 


m 4 


v 31 (m,/j) 


73i(m,/x) 


m 4 


to 4 


vw{m,n) 


733 (W, £t) 


m 4 


to 4 


U32(m,/i) 


732 (m,//) 


m 4 


to 4 



vi{m,n) +v 2 {m,n) +v 3 i(m,fi) +w 33 (m,^) + i; 32 (to,^), 
= -0.0506606 (2.log 



- 2.95381J , 

0.0000200507(.g^) 2 (183.83- 173.783 log (— ) - 48. U,-r 



Am/ 



M / 



0.0000228551 (^) 4 m 4 log 3 
0.0003291 14(<£) 4 ™ 4 log 3 , 
- {g Y ) 4 m 4 10- 3 (0.03291141og 5 ( - 



0.105904 log 4 ( — ) + 



rn 



0.0165851 log 3 (^j +0.441159 log 2 (^j - 0.74347 log (^j +0.388237). 



(86) 
(87) 

(88) 

(89) 

(90) 

(91) 
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Let us now consider that the renormalization point for /i is chosen at the same value of the fermion mass m/, 
under consideration, that is Log( 7 -^-) = 0. Also we will define new scaled scalar field $ and interaction parameter g 
by mean of 

$ = aip, (92) 
g y = am = gexp($), (93) 
g = am f , (94) 

Then, the evaluated total contribution to the effective potential for the Dilaton v(m, fi) can be expressed as a function 
w($, g) as follows 

v^g) _ vim^^ -0.0000329114e 8 V ( I )5 +0.000105904e 8 V<I )4 

raj raj, 

+0.000289673e 8 V$ 3 + e 4 * (-0.000441 159e 4< V - 0.000962436e 2 V ) $ 2 + 
e 4 * (0.00074347e 4 V + 0.00348448e 2 V - 0.101321) $ + 

e 4 * (-0.000388237e 4 V - 0.00368594e 2 V + 0.149642) . (95) 
Let us also define now the functions U5. U4 and U3 in the following form 
u 5 (®,g) _ v(<f>,g) 



4 4 
m f m, 



(96) 



7 "7 



4 

m 4 



0.000105904e 8 V* 4 + 0.000289673e 8 V^ 3 



e 4 * (-0.000441159e 4 V - 0.000962436e 2 V) $ 2 + 

e 4 * (0.00074347e 4 V + 0.00348448e 2 V - 0.101321) $ + 

e 4 * (-0.000388237e 4 V - 0.00368594e 2 V + 0.149642) , (97) 
" 3($ .'- 9) = + 0.000289673e 8 V* 3 + 



m 4 



e 4 * (-0.000441159e 4 V - 0.000962436e 2 V ) $ 2 + 

e 4 * (0.00074347e 4 V + 0.00348448e 2 V - 0.101321) $ + 

e 4 * (-0.000388237e 4 V - 0.00368594e 2 V + 0.149642) . (98) 

Note that U5 coincides v and is of order five in the powers of The function 1*4, U3 are defined as retaining only 
all the terms up to order $ 4 and <I> 3 respectively of the original function u 5 . Therefore, these functions basically 

correspond to the expansion of order five, four and three in powers of log (jjf^j ■ They are defined in order to study 

the influence of increasing the order of the perturbative expansion in powers of log (jfj ■ 

To evidence the dependence on $ and g of the three functions (after divided by the common factor rnA), they are 
plotted in figure [3] . The range of values of g — mja was chosen (0, 1) as suggested by the fact that a is of the 
order of the Planck length and thus the physical values of the considered fermion mass are expected to determine 
g to be smaller than one. The plot of U5 shows that there is a threshold value of g, below which the potential 
shows minima tending to stabilize the vacuum mean value of the Dilaton field. This behavior is also shown by the 
approximated potentials M4 and M3, a fact that indicates that after disregarding the higher quintic and quartic terms 
in the expansion in Log(^), the existence of Dilaton stabilizing minima is not affected. 

When considering the full evaluated potential curve M5 , illustrated at the top plot of figure G2 it can be observed 
that after lowering the g value below a critical threshold, the minimum as a function of $ stops to exist at a critical 
value g m in- However, in the case of 1*4 and 113 the minimum exists for arbitrary values of g— > 0. That is, when the 
potential approximations is bounded from below, the potential shows stabilizing minima at any small value of g close 
to zero. The field value at the minima grow when the coupling tends to vanish. It can be noted, that the non bounded 
from below character of the approximated potential calculated here is determined by the fact that the quintic power of 
$ correction turns to be negative. However, the physical system under consideration is one in which the total effective 
potential can be expected to show an exact bounded from below character. Thus, the next corrections are expected 
to exhibit a bounded from below behavior. In accordance with this expectation, in studying the g dependence at 
small values, we will employ the bounded from below approximated potential function U4, assuming that it represents 
a reasonably good approximation of the exact potential. 
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FIG. 3: The three figures show, from top to bottom, the potentials 115, 114 and M3 dependence on the field $ and the coupling 
g, respectively. The potential scale is chosen for a high magnification range (the minima of the surface at fixed g values are 
very far below the plotted range) in order to evidence the presence of a threshold for the appearing of the minima when the 
value of g decreases below g = 1. Note that for $ smaller than some units and not to small values of g, the three plotted 
graphs are similar, indicating that the elimination of the highest fifth, and also the next to highest fourth, powers of the field 
(or, of the logarithm in the original expansion) in defining 114 and 113 respectively, are not affecting the results in the mentioned 
region. The circumstance that the exact evaluated contribution has a negative leading term of order five (which makes the 
result unbounded from below) explains that for the plot of M5 the minima disappear for sufficiently small values of g. However, 
the fact that exact potential should be expected to be bounded from below, we consider that supports our assumption about 
employing the bounded from below approximations of the potential 114 in evaluating the Dilaton properties at the small values 
of g defined by the GUT and m tov mass scales. 



A. Dilaton field and mass for rrif at the GUT scale 



Let us consider now that the highest fermion mass m/ is given by the GUT mass scale 



m f = m GUT = 5.06773 x lO^cm -1 

= 10 16 GeV, (99) 

which produces for the coupling g the value 

3 

g = m f a = --K m GUT = -0.0030789542773. 

The potential u 4 as function of the field $ for this particular value of g is shown in figure 01 The minimum of the 
curve determines an estimate for the vacuum value of the Dilaton field given by 

<5> vac = 5.8576156 = a <p vac , (100) 

tfvac = -^5.8576156-. (101) 

3 K 
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FIG. 4: The effective potential ua defined by Eq. [97] as a function of the Dilaton field <£>. The fermion mass was fixed to 
correspond to the GUT mass itigut and the renormalization scale p, was chosen to coincide with this mass. The minimum of 
the potential is near the value $ = 5.7, which indicates that the field is bound to a high value near the Planck scale. 

This result indicates that the vacuum mean value of the Dilaton field, after assuming that the fermion mass is in 
the GUT scale, becomes stabilized in the scale of the Planck mass. 

Let us consider now the mass of the field excitation. Its value is determined by the second derivative of the potential 
curve taken at the minimum, which is given by 



it 4 ($, -0.0030789542773) 



1.28179 x 10 n m) 



In order to estimate the Dilaton mass let us consider the linearized equation of motion for the mean field 



1 d 2 
(^2<9 2 + ^« 4 ($, -0.00307895) 



)$ = 0, 



(102) 



(103) 



in which the factor multiplying the D'Alembertian appears due to the previously done change of field variable 
$ = a ip. 

The above wave equation leads to the dispersion relation for the Dilaton modes 



\p 2 + 4ztM®, -0.00307895) 
or a<tr 



( m GUT) 



.) =0 > 



(104) 



which for the case of the particle at rest p = (m^ lGUT \o, 0, 0) determines for the Dilaton the mass estimate 



ArnauT) 



d 2 

—uA§, -0.00307895) 
d<& 2 



m GUT I a 



5.58626 x 10 d 



(105) 



Therefore, the predicted order of the mass for the Dilaton also lays at extremely high values which make this 
field mode undetectable in a direct way. 



B. Dilaton mean value and mass for m/ at the top quark mass scale 

It is also of interest to take as mf the highest currently known fermion mass: that is, the top quark one 

.7164 x 10 15 cm" 1 . 



mtop = 172.0 ±0.9 GeV 

Then, the coupling g in this case takes the small value 

3 



g = nif (X 



-5.32659 x 10~ 17 . 



(106) 



(107) 
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FIG. 5: The effective potential U4 plotted as a function of the Dilaton field <E>. In this case the coupling was defined by a 
fermion mass correspond to the top quark one mt op and the renormalization scale fi was also chosen to coincide with this value. 
The minimum of the potential is now near the value $ = 36.7765, which indicates that the field is again staying at a high value. 



Figure [5] shows the dependence of the potential 114 as a function of the field $ at the above value of the coupling 
g. The minimum of the curve in this case gives for the mean Dilaton field at the vacuum 

$(,™^) = 36.3020096= a <p^f\ (108) 

V><"£°*) = -^36.3020096-. (109) 
3 K 

This result predicts that, assuming that the maximal fermion mass in Nature is given by the top quark one, which 
means a lower bound for the physical masses, the vacuum field of the Dilaton, again becomes stabilized in a scale, 
which although not being so high, is yet close to the Planck mass. 

In this case the dispersion relation for the Dilaton modes takes the form 



2 



{-\p 2 + ^o«4(*, -5.32659 x 10~ 17 ) 
or d<P z 



) = 0. (110) 



But, after evaluating for the second derivative of the potential at the minimum to be 

d 2 



-^uJ$, -5.32659 x 10~ 17> 



= 6.86404 x 10 64 , (111) 



and fixing again the rest frame momentum p — (m^J ltop \ 0, 0, 0) estimates for the Dilaton mass the value 



/ d 2 

m D ~~" = W ^2"4(*, -5.32659 x 10- 17 ) 



, X m top I a 

*=*;7c top 



= 7.07209 x 10^ y cm" 1 . (112) 

Henceforth, also in this case the predicted mass for the Dilaton turns to be a high value being now close to the 
GUT scale. Thus, it can be expected that for a maximal fermion mass in Nature ranging between the lower bound 
mtop and the GUT scale one, the Dilaton gets stabilized at a large field value as required by string phenomenology. 
In addition the resulting values of its mass, for the same range of rrif masses, is also out of the current observability 
range of particle detectors. 



V. CONCLUSIONS 



The predictions for the Dilaton stabilization problem determined the existence of massive fermion matter had been 
further investigated here. The fermion field mass values are considered in two cases: the top quark mass representing 
the lower bound of all existing but yet unknown fermion masses in Nature, and the energy scale of the grand unification 
theories of order 10 16 GeV. In both situations, the results indicate that the Dilaton mean field becomes stabilized at 
the very high values required by its role in allowing gravity to have its observed properties. Then, the same existence 
of matter seems to be a possible source of the dynamical fixation of the Dilaton field at the high values, required by 
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String Theory to imply the observable Einstein theory of gravity. Furthermore, the evaluations indicate that the 
Dilaton field is also found to be strongly bound around its mean value, by showing a large mass being close to the 
GUT or Planck scales. Therefore, the work identifies a possible explanation for the lack of observable consequences of 
the Dilaton scalar field in Nature. The discussion included contributions to the effective potential up to 3-loops. This 
allows to consider the influence of the inclusion of different leading perturbative correction on the main conclusions. 
After, disregarding in the evaluated potential: a) the highest order term (quintic) in the expansion in powers of 
Log(-) (which determined the unbounded from below structure of the potential at large $ values ) or b) the two 
highest orders (the quintic and the quartic ones), the obtained modified potentials are both bounded from below 
at high field values. This procedure allows that minima as functions of $ exist for arbitrarily small values of the 
coupling g. This allows to evaluate the small coupling values associated to the GUT and top quark masses. The fact 
that the Yukawa theory under consideration should exhibit a bounded from below potential, then supports this here 
adopted procedure for estimating the vacuum mean values and mass of the Dilaton field. However, further higher 
loop evaluations are convenient to define more precise estimated values of the Dilaton vacuum field and mass and also 
for checking that they do not affect the picture. The validity of the employed Yukawa like approximation in which 
the power expansion in Dilaton field of exponential factor of the Fermion mass is limited to the linear term in the 
field, is argued when the fermion mass is very much smaller than Plank mass. 
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Appendix A 

In this Appendix we discuss two points which improve the argues in the work: a) Adding a discussion about the 
validity of the Yukawa approximation; and b) To clarify the exposition of the results which indicate the stabilization 
of the vev and the mass generation. 

Below we address the two main critical remarks we identified in the Report: 



1. Yukawa approximation 

Let us study the validity of only retaining the linear term in the expansion of the exponential of the Dilaton field 
multiplying the fermion mass in the action. For this purpose let us consider the mean value of the square of the 
radiation field 4> r times the Dilaton coupling a in the just considered Yukawa model. The square root of this 
quantity gives a measure of the amount of the quantum fluctuations of a 4> r (x). A resulting small value of this 
quantity justifies the approximation done in the work. In terms of the Dilaton Green function D this quantity can be 
written as follows: 



a?iD*fo(»)|0>=afo / J^rMP) 



2 

a (R) 



dp d 
dp d 



(2-K) d i -p 2 + U{p,m,d) 



2 f dp d e 1 ni „, 

- a ^J (2ny P i+ii(p e ,m ) dy [116) 

where LT is its selfenergy which will be evaluated in the one loop approximation, by excluding the divergent pole parts. 
This corresponds to using the MS substraction by also employing the Feynman expansion in terms of the renormalized 
fields incorporating counterterms. 

Let us now consider that the bare Dilaton theory being quantized, is valid in the Planck scale where its action 
is assumed to define a low energy physics for a string theory. Then, after denoting with the B subindex the bare 
quantities, the relations between the bare and renormalized magnitudes are 

a (B) = ~k = Z%af R) , (114) 

d>\ B) (x) = zlct>l R) (x). (115) 

We will also assume that the product a 2 R ^ 4>(r) i x ) > after the study of the renormalization properties of the theory, 
can be chosen as a renormalization group invariant, which determines for the renormalization constants the relation 

z*- 1 

Zj a — — — . 

Zl 
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Therefore, the squared deviation of the field can be written in the form 



«Cfl)^(B)(«)|0) = -^ i <0| ^(z) |0) 

= *ZL{J*L i (us) 

Now, in order to get an estimate of this quantity, let us consider the approximation in which the renormalized Yukawa 
coupling is small and also assume that the finite part of the one loop selfenergy is evaluated at zero momentum. As 
it will be seen in the comment of the next point, this quantity determines the first approximation for the Dilaton 
mass as given by the square root of the value of p\ fixing the pole of the propagator. As it was mentioned, in 
the discussion of the next point, assuming that the coupling is small implies that the selfenergy evaluated at zero 
momentum approximately determines the pole mass of the Dilaton. Thus in this small coupling limit at least the above 
described approximation can give a reasonable estimate of the field fluctuations. In addition, as it was supposed in the 
work, we will fix the scale parameter \l as coinciding with the fermion mass, that is fJ, = m. Therefore, approximating 
IT 1 ) (p e , m, d) by il*- 1 ' (0, m, d) and evaluating the remaining simple momentum integral gives 

"(B) 



In order to proceed, let us evaluate the wavefunction renormalization constant Z$ to further transform the previous 
formula. The fermion one loop contribution without the substractions is given by the expression 

dq d 1 1 



U(p,m,d) = -gl J j0sDfr[ 



(2n) D i m — j^p^ m — 7^(p — q)^ 

to which should be added the wavefunction and mass first order counterterms in order to evaluate IT 1 ) in formula 
(|117p . Calculating the Dirac traces and performing the Wick rotation allows to explicitly determine II as follows 

dq d e 1 



IL(p e ,m,D) = -Ag 2 y m d - 2 J 



2 2 
9y P, 



(2ir) d 1 + & 
dq d 1 1 



v I e J (2ir) d l + q 2 e l + ( Pe - q e ) 
Til - -) 



2 



(2tt)» 

St Pl'^W 1 f d < 1 + 4^(! " ( 118 ) 

(47TJ2 Jo m 

The Feynman parametric integral can be explicitly performed by using the Wolfram Mathematica program in the 
form 



.F(4>e)= [ dx(l + ^x(l-x)y 
Tn Jq m 

2-( 1+£ )(g + V¥+4)(i 



2^(1-6,6,2-6,1(1+^=1==))- 
2-(l+=)(g_ v ^ + I)(l + 



^(1-6,6,2-6 ,1(1+—*)), (119) 



Re[^±i > 1] or Re[^±i < 1] or * real. (120) 
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The exact expression for the finite part of the one selfenergy can be obtained by summing the counterterms. This 
permits to define the wave function renormalization constant and the coefficient of the mass counterterm as the pole 
part in e of the fermion loop contribution. The wavefunction renormalization constant is given as 

This formula allows to find the limit d = 4 — 2e — >4 limit of the mean of the squared deviation of the field in 
(fTTT)) as follows 



^(7r) 2 -r( e -i) 



nW(o,m,^ <t;^ i; 



2 ix "*>"7 (2tt 

= a lim 



i=?6 i (<#) 2 1 

= 167r 2 (l-7)a 2 m 2 

= i6(_)V(l- 7 ) k 2 ™ 2 , (121) 
7 = 0.57721. (122) 

This result indicates that the mean squared deviation of the argument of the exponential of the Dilaton radiation 
field is of the order of the square of ratio between the mass of the fermion and the Planck mass. Therefore, for the 
values considered in the work: the top quark mass and the GUT unification scale, the exponential interaction term 
in the Dilaton radiation field should be well approximated by the linear term defining the Yukawa model. 

2. Approximations in determining Dilaton stabilization and mass 

In this subsection, we describe the main criteria of stabilization of the vev and the approximations done in evaluating 
the mass in the paper. Let us consider the definition of the effective action of a scalar field 

T[cj>{x)] = -MZ\j]- J dx j(x)ct>(x), 

and its expansion around an homogeneous vev (f> in a functional power series in spacetime dependent fluctuation f{x), 

T[(j) + ip(x)} = T[<j)} + J dxj(x)<p(x) - ^ J J dxdy(p(x)D^ 1 (x - y)ip{y) + ... 

= -v^v[cb] + J d P f(p)<p(p) -^[[ ^(p)d- 1 (pMp) + ... 

= V^(V[4>} + I dpf(p) V (p)- f f ^(p)D- 1 (pMp)...) 7 (123) 



where due to the homogeneity of the vev T[4>] is proportional to the four dimensional volume V^ D \ The negative 
sign in the quadratic in the field term of the expansion comes because the second functional derivative of the effective 
action respect to the field is the negative of the inverse of the propagator. The quantity V appearing defines the so 



called effective potential V[4>\ — — . The potential V has the interpretation of the energy density of the quantum 
system at the value of the also homogeneous external source that sustains it. The stable states of the vev fields are 
given by the minima of V [20I l2lj ] . This property implies that the results in the work indicate the stabilization of 
the Dilaton vev due to the quantum corrections associated to massive matter fields. 

Lets us consider now the question of the approximation in which the results determine an estimation of the mass 
for the oscillation around the mean field. The discussion in the previous subsection partially helps to clarify this point. 
The masses of the particles are defined by the squared momenta making equal to zero the inverse propagator 

D-\p) = -{p 2 -n (1) (p,m,4)) = 0. 

where, since II^ 1 ) is finite, d = 4 has been substituted. The approximation which is employed to estimate the Dilaton 
mass in the work is obtained by expanding the proper mass (the squared momentum value) which makes vanish the 
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inverse propagator, in powers of the coupling constant gy as follows 



P 2 (9y)= E M) 



m=0,l,2... 



After substituting this series in the inverse propagator, and noting that the inverse propagator for the scalar field is a 
function of the momentum only through its squared values and that the coupling expansion of the one loop selfenergy 
is already of order (gy) 2 , it follows 



Therefore, in the first approximation, the proper mass is given by the squared root of the zero momentum component 
of the selfenergy. Further, the selfenergy at zero momentum coincides with the squared root for the second derivative 
of the effective potential, and this was the criterion employed to evaluate the Dilaton mass in the work. This property, 
can be seen after expressing the effective action expanded around an assumed minimum and homogeneous field value 
as follows 



-Po = °. 

(-p?(.9?-) 2 + n(o,™,4)) = o, 



m 2 =p 2 =n< 1 >(0,m,4) + O(G&) 4 ). 



(124) 




+ ... 



V.(4)y M __^_ [£) -l ( p)] p=0 + ... 

v^(v[^ + ^[n( P )] p=0 + ..) 
v^ivw + ^-^vw + ...,) 



(125) 




P(p)] P =o. 



(126) 



